Efficient Spin Squeezing with Optimized Pulse Sequences 
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Spin squeezed states are a class of entangled states of spins that have practical applications to 
precision measurements. In recent years spin squeezing with one-axis twisting (OAT) has been 
demonstrated experimentally with spinor BECs with more than 10^ atoms. Although the noise is 
below the standard quantum limit, the OAT scheme cannot reduce the noise down to the ultimate 
Heisenberg limit. Here we propose an experimentally feasible scheme based on optimized quantum 
control to greatly enhance the performance of OAT to approach the Heisenberg limit, requiring 
only an OAT Hamiltonian and the use of several coherent driving pulses. The scheme is robust 
against technical noise and can be readily implemented for spinor BECs or trapped ions with current 
technology. 
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Spin squeezed states [ij have attracted a lot of interest 
due to both its role in the fundamental study of many- 
particle entanglement and its practical application topre- 
cision measurements with Ramsey interferometers [2|-l6j- 
In recent years, much progress has been made on the ex- 
perimental squee zing of a large number (10'^ ~ 10^) of ul- 
tracold atoms |7Hll| . Many of these experiments follow 
the so-called one-axis twisting (OAT) scheme, which is 
known to reduce the noise/signal ratio from the classical 
case by a amount that scales as A^~^/'^with the particle 
number N [1\. This reduction is not optimal yet and still 
above the so-called Heisenberg limit which scales as N~^. 
There have been several theoretical proposals to enhance 
the OAT (l3 . [l3fl . For example, one of the approaches 
[13| involves inducing a better squeezing Hamiltonian, 
the so called two-axis twisting (TAT) Hamiltonian, with 
Raman assisted coupling for trapped spinor BECs. This 
is a hardware level engineering, requiring modification 
of a particular experimental setup and does not apply 
to other physical systems. Another approach [l^l em- 
ploys a digital quantum simulation technique to convert 
an OAT Hamiltonian to an effective TAT Hamiltonian by 
stroboscopically applying a large number of pulses. This 
software level solution is universal but sensitive to the 
accumulation of control errors. None of these propos- 
als have been experimentally tested yet due to various 
difficulties. 

Inspired by the idea of optimized quantum control, we 
propose an experimentally feasible scheme to greatly im- 
prove the performance of OAT, requiring only two or 
three additional coherent driving pulses to carry out col- 
lective spin rotations, which is a routine technique with 
the current technology. The scheme is shown to be ro- 
bust to noise and imperfection in control pulses. Using 
this scheme, it is possible to generate more spin squeezing 
and detect a significantly larger entanglement depth for 
the many-particle atomic ensemble [5|. This new scheme 
enhances the OAT squeezing on the software level and 
therefore can be applied to any physical system that is 
endowed with these operations. The idea of optimized 
squeezing may also be easily transferred to cases where 



the interaction term deviates from the OAT Hamiltonian. 

We consider the general scenario of one-axis twisting 
independent of the underlying physical system with the 
Hamiltonian H = xSl {Sz = Y.f K) (setting h = 1). 
The system starts from a collective spin coherent state 
polarized along x-axis. As time goes on the initially 
homogenous spin fluctuation gets distorted and redis- 
tributed among different directions and the direction 
along which spin fluctuation gets suppressed gradually 
changes over time. The squeezing is measured by the 
parameter ^^, defined as ^^ = N (S?^) / \{Sx)\ , where 
n is the direction along which spin fluctuation is mini- 
mized. The decreasing rate of ^^ slows down with time, 
and after the optimal squeezing point, ^^ increases again. 
Aside from the initial state, which is rotationally sym- 
metric about a;- axis, all the subsequent states breaks this 
symmetry and picks out a special direction, i.e. the 
direction along which fluctuation is minimized. It is 
well known that the two-axis twisting (TAT) Hamilto- 
nian H2 = X2 (S'^ — Sy) can produce better squeezing 
[1'], which, after doing the Trotter decomposition with 
an infinitesimal time interval, could be seen as switching 
the squeezing axis back and forth very fast between two 
orthogonal directions [IJI- To avoid the noise accumu- 
lation from a large number of switching pulses inherent 
in the Trotter expansion scheme, we take an alternative 
approach based on optimization of a few control pulses 
to maximize the squeezing of the final state. We consider 
an n-step squeezing protocol (where n is typically 2 or 
3 for a practical scheme) defined as follows: at step j 
{j — 1, 2, ..., n), we first apply an instantaneous collec- 
tive spin rotation around x-axis, U{ai) — exp{—i SxOn), 
and then let the state evolve under the OAT Hamiltonian 
H = x^l for a duration T^. Effectively, we squeeze the 
state along a different axis lying in the y — z plane in each 
step, so the effective evolution operator can be written 



U{e,,T,) = \{exp{-ixSlT{), 
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Figure 1: The squeezing ~log{^^) as a function of the control 
parameters O2 and T2 for a typical value of Ti, calculated with 
N=2000 spin-1/2 particles. See Eq. [T]and text for definition 
of 6i and Ti. The cross symbol marks the point of optimal 
squeezing. The horizontal line ^2 = corresponds to the case 
of the OAT scheme. 



where Sg = cosOjSz + sinOjSy and the factors are ar- 
ranged from right to left with increase of j. Since the 
initial state is assumed to be polarized along x-direction, 
which is symmetric around z-axis, 9i is irrelevant and can 
be chosen to be (so no control pulse is needed for step 
1). Therefore, for an n-step squeezing protocol, there are 
{2n — 1) tunable parameters: Ti and 9i (excluding 61). 
The final squeezing parameter is thus a multi-variable 
function ^'^{Ti, 9i). Our purpose is to find the best avail- 
able squeezing ^^(T^, 9i) with a minimum number n of 
the time steps. 

In the case of n = 2 or 3, the landscape of £,'^(Ti, 9i) in 
the parameter space is quite simple and well behaved. 
Take the n = 2 case as an example. For a typical 
value of Ti smaller than the optimal OAT squeezing time, 
— /og(^^) as a function of 6*2 and T2 is shown in Fig[TJ The 
optimal squeezing point marked by the cross lies way off 
the OAT trajectory, the horizontal line with 6*2 = 0. For 
the n = 3 case, with 02 and T2 fixed near the optimal 
values of the n = 2 case, —log{£,'^) as a function of ^3 and 
T3 shows a similar landscape. These solutions already 
exceed that of the OAT scheme by a large margin. The 
results indicate that the optimization technique with n 
as small as 2 or 3 suffices to significantly improve over 
the OAT scheme. 

Next, we investigate performance of the optimized 
squeezing scheme, focusing on the scaling of the squeez- 
ing ^^(Ti, 9i) as a function of the total particle number 
iV. For a given set of parameters, we can numerically 
calculate the evolution operator in EqU] by exactly diag- 
onalizing the effective Hamiltonians Sg, and then obtain 
the squeezing parameter ^^. We randomly sample from 
the parameter space for a large number of times, use 
these random samples as initial guesses to start uncon- 
strained local optimization of the squeezing parameter, 
and pick the best one as our solution. Repeating this 
procedure for every system size TV is extremely resource 
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Figure 2: Scaling of the squeezing parameter ^^ with the num- 
ber of qubits. Curves from top to bottom are for one-axis 
twisting (OAT) , two-step optimized squeezing, two-axis twist- 
ing (TAT) and three-step optimized squeezing. Inset shows 
the same curves in log-log scale. 



intensive especially when N gets as large as 10^. Taking 
advantage of the fact that adding several more to 10'^ 
particles should not change the solution much, we can 
feed the previously found non-local optimal solution as 
an initial guess to the local optimizer of a larger system 
and obtain a near optimal solution quickly. In this way 
we managed to obtain (near) optimal solutions for sys- 
tems all the way up to iV = 10^ particles, with only a 
cost of classical computing time on the order of tens of 
hours on a typical multi-core computer. As shown in Fig 
[21 with n — 2, the squeezing parameter ^^ gets reduced 
by a significant amount already compared with the OAT 
scheme, and with n — 3, £,^ decreases further. The scaling 
of ^^ with the number of particles shows a clear power 
law ^2 ^ l/Nl^. A simple OAT scheme gives P = 2/3 
and the TAT scheme gives /3 = 1 [1|. The Heisenberg 
limit of noise gives a bound /3 < 1 for the scaling, and 
this bound is saturated by the TAT scheme. Remarkably 
we observe that the optimized n = 2, 3 protocols can give 
/3 = 0.92 and 0.98, respectively, very close to the ultimate 
Heisenberg limit. Moreover, the n — S optimized scheme 
has a smaller multiplicative constant compared with the 
TAT scheme, so in the realistic range of particle number 
N < 10^, it actually outperforms the TAT scheme. This 
shows that a moderate alternation of the OAT scheme 
through optimization can significantly increase the spin 
squeezing. 

We have demonstrated a significant improvement over 
the conventional OAT by applying very few optimized 
control pulses. A cost of the proposed scheme is that it 
takes longer evolution time to achieve the optimal squeez- 
ing. A typical evolution of ^^ with time t is shown in 
Fig. [3l We notice that in general the {i + l)-th squeez- 
ing step takes longer time than the i-th step. Since the 
time cost in the first step is on the order of the optimal 
OAT duration, the overall duration of the new protocol 
is usually longer than that of the OAT scheme. An exces- 
sively long duration would be an obstacle in systems with 
short coherence time. The two relevant time scales here 
are the coherence time r and the inverse of interaction 
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Figure 3: Evolution of tiie squeezing parameter ^'^ with time, 
calculated with N=2000 spin-1/2 particles. The dash-dot line 
is for one-axis twisting (OAT), the dash line for the two-step 
optimized squeezing scheme, and the solid line for the three- 
step optimized squeezing. 
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Figure 4: Constrained optimization of ^^ with the total time 
duration as a cost function. We take 1/x as the time unit. 
Achievable squeezing ^^ as a function of the total duration 
is shown, together with one-axis twisting (OAT), calculated 
with iV = 2000 spin-1/2 particles. OPT-2 (3) stands for opti- 
mized squeezing sequence with n — 2(3) segments. Horizontal 
and vertical dashed lines are guides to the eye. 



strength l/x- The time cost of the new scheme is around 
0.01/x ~ 0.1/x- If ''■ ^ 0.1/x the new scheme can be 
implemented without compromise. On the other hand, if 
that is not the case, decoherence effect would play a role 
and our unconstrained optimization no longer yields the 
best result. However, we can work around this problem 
by performing an optimization with the total duration 
added as a cost function and get a compromised optimal 
pulse sequence. By tuning the weight of the cost func- 
tion we could obtain a continuous series of compromised 
optimal solutions as shown in Fig. 21 These solutions of 
two-step and three-step schemes form two line segments, 
continuously connecting the optimal OAT squeezing pro- 
tocol to that of the unconstrained optima, offering a trade 
off between the protocol duration and the squeezing mag- 
nitude. For each real experimental setup, one could cor- 
respondingly pick up the best point in accordance with 
the coherence time of the system. How much one can gain 
over the OAT scheme depends on how long the coherence 
time can reach. 

Next we test noise resistance of the proposed scheme. 





Figure 5: Optimized squeezing in the presence of control 
noise. We use the three-step optimization scheme as an exam- 
ple and assume all the five control parameters in this scheme 
have the same magnitude of relative errors as specified in this 
figure. The dash line is for the ideal case with no error in 
the control parameters, the solid line denotes the average of 
many random trajectories (about 50 random trials) and the 
shaded area marks the range of those trajectories. In the left 
panel, the shaded region is too small to be distinguished from 
the ideal case. 



There are only 3(5) control parameters in the n = 2(3) 
scheme, making the accumulation of control noise negli- 
gible. We have done numerical simulation of our scheme 
adding random pulse area/timing noise and confirmed 
the robustness of the squeezing parameter ^^ as shown 
in Fig. [5l This contrasts to the proposals (l2 . [lJ| re- 
quiring a large number of coherent rotation pulses where 
control errors accumulate and significantly degrade the 
performance. Thus our proposed scheme offers a useful 
alternative to the previous works. Another practical is- 
sue related to control noise is the uncertainty in number 
of particles in a real experiment. Our pulse scheme de- 
pends on the number of particles TV while in experiments 
such as ultracold gas we do not typically know the num- 
ber N exactly. Fortunately we notice that the control 
parameters vary slowly with N and an uncertainty in N 
is equivalent to a small extra noise in the control param- 
eters, to which ^^ is not so sensitive as we have shown in 

Fig. [g 

Finally we discuss possible physical realizations of the 
scheme proposed here. The scheme only requires two 
ingredients, the nonlinear collective spin interaction S'^ 
and the ability to rotate the collective spin around an or- 
thogonal axis, say x. Several experimental systems meet 
these requirements, e.g., trapped ions and spinor BECs. 
In trapped ion systems, depending on the ion species, one 
can use bichromatic lasers or two pairs of Raman laser 
beams (the Molmer-Sorensen scheme) to induce the S^ 
or S'^ type of interaction. The strength of this inter- 
action X can reach kHz scale, giving l/x '^ ins. The 
coherence time usually exceeds l/x and our scheme can 
apply without compromise. Collective spin rotation can 
be simply done by shining laser on all the ions driving the 
corresponding single-qubit T^/y or rotation. The rotation 
pulses have durations much shorter than 1/x- While lin- 
ear Paul traps [13| can now coherently control only about 
a dozen of ions, too few for the purpose of spin squeez- 
ing, planar Penning traps can manipulate more than 200 
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Figure 6: The entanglement depth achievable with different 
approaches for 200 spin- 1/2 particles. The solid lines from 
top to bottom correspond respectively to the OAT scheme, 
the two-step optimized squeezing, the TAT, and the three- 
step optimized squeezing. The dashed lines from top to bot- 
tom correspond to the optimal squeezing for 50, 100, and 
200 particles respectively. Lying below the curve of optimal 
squeezing for n particles is a certificate of genuine n-particle 
entanglement. 



ions [161 . For the purpose of precision measurement, 200 
ions may seem less impressive than 10^ particles, but we 
show that using our scheme we can create genuine multi- 
particle entangled states with a significantly larger entan- 
glement depth. The entanglement depth, defined in [5J, 
is a way to measure how many particles within the whole 
sample have been prepared in a genuine entangled state. 
Our result is shown in Fig. [6l In this figure, a point lying 
below the optimal squeezing curve of n particles corre- 
spond to a state that contains genuine n-particle entan- 
glement. Our scheme produces states that lie below the 



OAT states in a large range of (Sx) values, which means 
that experimentally one can achieve a significantly larger 
entanglement depth by this optimization technique. 

Another class of physical system is a spinor Bose- 
Einstein condensate of atoms with two chosen internal 
states mimicking spin- 1/2 particles [^,l3|. The desired S'^ 
interaction is induced by spin-dependent s-wave scatter- 
ing as proposed in [4J . Coherent laser pulses illuminating 
the whole condensate can implement spin rotations sim- 
ilar to the trapped ion case. However, the strength of S'^ 
interaction is much smaller compared with the trapped 
ion case, x = 0-3 ^ 0.5 Hz as reported in [3, Q. The 
coherence time for the spinor BEG is also shorter. Hence 
we typically need to apply the compromised scheme, us- 
ing the actual coherence time and interaction strength of 
the system as input parameters. 

In summary, we have proposed a new method based 
on optimization to significantly enhance spin squeezing 
using the one axis twisting Hamiltonian. To achieve sig- 
nificant improvement in spin squeezing, we need to apply 
only one or two global rotation pulses at an appropriate 
evolution time and with optimized rotation angles. Us- 
ing two pulses, the final squeezing is very close to the 
Heisenberg limit already. As we use a very small number 
of control pulses, the scheme is immune to accumulation 
of control errors and can be readily applied in experimen- 
tal systems without significant modification of the setup. 
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